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We examine the conditions needed to accomplish stimulated Raman adiabatic passage (STIRAP) 
when the three levels (g, e and /) are degenerate, with arbitrary couplings contributing to the 
pump- pulse interaction (g - e) and to the Stokes-pulse interaction (e-f) . We show that in general a 
sufficient condition for complete population removal from the g set of degenerate states for arbitrary, 
pure or mixed, initial state is that the degeneracies should not decrease along the sequence g, e and 
/. We show that when this condition holds it is possible to achieve the degenerate counterpart of 
conventional STIRAP, whereby adiabatic passage produces complete population transfer. Indeed, 
the system is equivalent to a set of independent three-state systems, in each of which a STIRAP 
procedure can be implemented. We describe a scheme of unitary transformations that produces 
this result. We also examine the cases when this degeneracy constraint does not hold, and show 
what can be accomplished in those cases. For example, for angular momentum states when the 
degeneracy of the g and / levels is less than that of the e level we show how a special choice for 
the pulse polarizations and phases can produce complete removal of population from the g set. Our 
scheme can be a powerful tool for coherent control in degenerate systems, because of its robustness 
when selective addressing of the states is not required or impossible. We illustrate the analysis with 
several analytically solvable examples, in which the degeneracies originate from angular momentum 
orientation, as expressed by magnetic sublevels. 

PACS numbers: 32.80.Qk,42.65.Dr,33.80.Be 



I. INTRODUCTION 

Techniques based on adiabatic passage provide very 
practical methods for producing nearly complete transfer 
of population between two quantum states using crafted 
laser pulses [1] . One popular example of such coherent 
adiabatic excitation, stimulated Raman adiabatic pas- 
sage (STIRAP) [2], provides a simple and robust tech- 
nique for transferring population between two nondegen- 
erate metastable levels, making use of two pulses, termed 
the pump pulse (linking the initially populated ground 
state ip g with excited state ip e ) and the Stokes pulse (link- 
ing excited state ipe with final state ipf of the three-state 
chain). When the pulses are properly timed (Stokes pre- 
ceding but overlapping the pump pulse) and two-photon 
resonance is maintained, then via adiabatic passage the 
population is transferred from initial to final state, with- 
out appreciable population in the excited state at any 
time. 

The operation of STIRAP can be understood by in- 
troducing instantaneous eigenstates of the time-varying 
Hamiltonian, the time-dependent adiabatic states with 
associated time-dependent eigenvalues (adiabatic ener- 
gies). One (and only one) of these states, $o(Oj is con- 
structed from only the initial and final state, with no 
component of the excited state. Because the excited state 
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generates fluorescence via spontaneous emission, such an 
adiabatic state will exhibit no such signal; it is termed 
a dark state. During the STIRAP process the state vec- 
tor ^f(t) remains aligned with the adiabatic state &o(t), 
while this state, in turn, changes composition from being 
aligned with ip g initially to being aligned with ipf after 
the Stokes-pump pulse sequence. 

Numerous extensions of the basic three-state STIRAP 
[2] have been considered [1, 3], including examples in 
which there occur magnetic sublevels and associated de- 
generacy. One possibility is that the atomic energy levels 
are coupled in such a way that each one is connected to 
at most two others. Population transfer in such multi- 
state chains has been studied by several authors [4-12]. 
In addition to straightforward population transfer, STI- 
RAP has been applied to the problem of manipulat- 
ing and creating coherent superpositions of two or more 
quantum states. Such superpositions are required for 
many contemporary applications including information 
processing and communication. The original STIRAP 
process has, for example, been utilized to create coherent 
superpositions in three- and four-level systems [13-19] 
and to prepare A-component maximally coherent super- 
position states [20] . There have been proposals to create 
TV-component coherent superpositions in such systems, 
where the final state space is degenerate [21, 22], at least 
in the rotating wave picture. This idea has been further 
developed to map wave-packets between vibrational po- 
tential surfaces in molecules [23, 24]. Finally, it has been 
shown for a specific degenerate system, having a single 
initial-, two degenerate intermediate-, and three degener- 
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ate final states coupled in the Raman configuration, that 
the STIRAP process can be extended to systems with 
degenerate intermediate and final levels [25]. 

Yet an open question has remained: what is the most 
general system of three degenerate levels, linked via Ra- 
man process, for which it is possible to transfer all pop- 
ulation from the ground-state manifold of degenerate 
states (the g set) to the final-state manifold (the / set) 
while minimizing population in the excited states (the 
e set), without first using optical pumping to prepare a 
single nondegenerate initial state? We here provide the 
answer to this question. 

We consider N g degenerate states of the g set, coupled 
by means of a pump-pulse to N e degenerate states of the 
e set, which in turn are linked by the Stokes pulse to Nf 
degenerate states of the / set. We will show that such 
a generalized STIRAP process is almost always possible 
if the succession of state-degeneracies is nondecreasing, 
i.e. N g < N e < Nf. When such conditions hold, then 
for arbitrary couplings among states (e.g. arbitrary el- 
liptical polarization of electric dipole radiation between 
magnetic sublevels) it is possible to obtain complete adi- 
abatic passage of all population from the states of the g 
set into some combination of states of the / set. 

We also examine the possibility of adiabatic passage 
when this restriction on degeneracies does not hold. We 
show that in this case in general only part of the popula- 
tion can be transferred to the / set. We point out that, 
in special but important cases, for an appropriate choice 
of the polarizations and phases of the coupling fields, a 
complete adiabatic population transfer can be obtained. 

Another motivation of this paper is the creation of co- 
herent superposition states in a degenerate system. The 
difficulty in such systems arises from the limited possibil- 
ity of addressing a single preselected state: addressing of 
a selected state is usually achieved by exploiting selection 
rules that the coupling field should satisfy. However, if 
we have e.g. two Zeeman multiplets a light field with a 
certain polarization will create several couplings between 
the magnetic sublevels of the multiplets. Our scheme of- 
fers a solution to this problem: we show that despite of 
the lack of selective addressing of the degenerate states, 
we have some control over the created coherent superpo- 
sition state in the / set. As we point out, and illustrate 
with specific examples, the level of control depends on 
the system under consideration. 

Our scheme is based on using a Morris-Shore (MS) 
transformation of the Stokes couplings or the pump cou- 
plings, thereby reducing this particular (generally com- 
plicated) linkage to a set of unlinked two-state systems 
and dark states [26, 27]. Underlying this technique is 
the fact that, as Morris and Shore [26] have shown, any 
system of linkages in which there occur only two detun- 
ings (i.e. the system has two sets of degenerate sublevels, 
termed here a and 6, forming sets of dimension N a and 
Nt), can be transformed, via suitable redefinition of basis 
states, to one involving a set of A< independent two-state 
systems, where N < = min{A a , Nb}, together with a set 



of uncoupled states that are unconnected to other states 
by the given couplings (one-state systems). If such an un- 
coupled state has no component from the e set we term 
it a dark state. We here extend that work to produce 
sets of unlinked three-state systems. 

The paper is organized as follows: In the next sec- 
tion we present a general model for degenerate, three- 
level systems and discuss its main properties. In Sec. Ill 
we derive a general condition for complete STIRAP-likc 
population transfer. In Sec. IV we derive analytic ex- 
pressions for the dark and bright states for important 
special choices of degeneracies. Then, in Sec. V, we de- 
termine the conditions needed for adiabatic evolution. 
We demonstrate our method through some specific ex- 
amples in Sec. VI. Finally, in Sec. VII, we summarize 
our results. 



II. THE DEGENERATE-SUBLEVEL MODEL 



A. The Hamiltonian 

As is customary when dealing with STIRAP or other 
three-level chains, we introduce an expansion of the state 
vector ^(t) that incorporates explicit phases taken from 
carrier frequencies of the pump and Stokes pulses, to p and 
los, respectively. In this rotating- wave picture, and with 
the customary neglect of counter-rotating terms [i.e. time 
variations (tOi + tOj)t] the rotating- wave approximation 
(RWA) Hamiltonian takes the block-matrix form 



H(t) 



p(t)P 
p(t)Pt hA s(t)S 
«(t)St o 



for use with the Schrodinger equation 



ih— C(t) = H(t)C(t) . 



(1) 



(2) 



Here the zeros denote null square or rectangular ma- 
trices of appropriate dimensions. The zero matrix in the 
bottom right corner indicates that the system is supposed 
to maintain two-photon resonance. All time dependence 
occurs in the two pulse amplitudes p(t) and s(t), each 
with unit maximum value. The N e x N e diagonal matrix 
fiA describes the detuning of the pump carrier frequency 
from the Bohr frequency of the g — e transition. The 
N g x N e matrix 2p(i)P/H consists of Rabi frequencies 
associated with the transitions between the g and e sets, 
HSlij (t) = 2p(t)Pij . The elements of the constant matrix 
P read 



P 



i = 1 • • • Ng 

3 = l---N e 



(3) 



where £^ is the peak amplitude of the pump-pulse elec- 
tric field and fiij is the dipole-transition moment. 

Similarly, the N e xNf matrix 2s(t)S/h consists of Rabi 
frequencies associated with the transitions between the e 
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FIG. 1: (Color Online) An example for the degenerate STI- 
RAP scheme: we have three Zeeman multiplets with J g — 2, 
J e = 3, Jf = 4. The couplings are those of a± polarized 
pulses. The pump and Stokes pulses are detuned from ex- 
act resonance with the excited-state by A, but they maintain 
two-photon resonance between states g and /. The system 
separates into two independent systems, indicated by solid 
and dashed lines. 

and / sets of states. The elements of the constant matrix 
S are 

&ij = 2^ ' I j — l . . . TV/ ' ^ 

where £^ is the peak amplitude of the Stokes electric 
field. 

The structure of the RWA Hamiltonian of Eq. (1) is 
similar to that of the conventional three-state STIRAP, 
in having all time dependence confined to two pulses p(t) 
and s(t), but instead of single ground, excited, and fi- 
nal states we have degenerate manifolds of sublevels, and 
hence we have matrices p(t)P, s(t)S, and A where con- 
ventional STIRAP would have scalar elements. To illus- 
trate these Fig. 1 shows the linkage patterns for the an- 
gular momentum sequence J = 2 <-> 3 <-> 4. To simplify 
the drawings we show the energies of successive manifolds 
as increasing, such as would occur with a ladder scheme; 
the connections are the same as with the usual lambda 
couplings, in which the final sublevels have energies be- 
low the excited state. 

Although we discuss situations in which the coupling 
matrices result from magnetic-sublevel degeneracy, all of 
our results apply quite generally, for any mathematical 
form of the dipole-moment matrices and consequently for 
any arbitrary structure of the constant matrices S and 
P. 



B. Dark states 

There exist N = N g + N e + Nf basis states for this 
system, and hence adiabatic states We can 

immediately apply the MS transformation [26], at each 
instant of time, by placing the g and / sets of states 
together into the MS a set, and taking the e set to be the 
MS b set. If the a set is larger than the b set, there will be 
N u = N a — Nb uncoupled states. None of these have any 



component from the e set, and so they are all dark states. 
The number of dark states is thus Njj = N„ + Nf — N e . 
In the conventional nondegenerate STIRAP [3] , for which 
N = 3, the MS transformation gives one dark state and 
one bright state; for the tripod system, for which N — 
4, there are two dark states [17, 18]. In the angular- 
momentum system of Fig. 1 there are Nd = 5 + 9 — 7 = 7 
dark states. 

For conventional nondegenerate STIRAP the composi- 
tion of the dark state changes with time, because the cou- 
pling matrices and the MS transformation change with 
time. However, it is possible to associate the (single) dark 
state initially with the nondegenerate ground state by ap- 
plying the pulses in the counterintuitive order, i.e. Stokes 
pulse preceding pump pulse. When there is degeneracy, 
it is necessary to establish that the entire population of 
any pure initial state in the g set is projected into the set 
of dark states and no population is left in bright states. 
This completeness of the dark states is at the heart of 
our question concerning the possibility of STIRAP with 
degeneracy. 



III. GENERAL CONDITION FOR COMPLETE 
POPULATION TRANSFER 

One of our basic questions is whether, for a given link- 
age pattern, it is possible to empty completely the g set 
for any arbitrary initial state, once we have fixed the 
pump and Stokes pulses. 

It is easy to see that one necessary condition for com- 
plete removal of population from the ground manifold 
is that there should not be more sublevels in this man- 
ifold than there are in the excited state, i.e. we require 

Ng<N e 

To prove this assertion we employ a MS transformation 
[26] on the pump transitions that connect ground and ex- 
cited states. This transformation introduces a new set of 
basis states in each of these manifolds, such that each 
sublevel from the g set couples to at most one sublevel 
from the e set. Were there are no Stokes couplings be- 
tween e and / states, the dynamics could be described 
as a set of independent two-state systems, together with 
some single states (uncoupled states) that are not af- 
fected by the pump radiation. Given such a revision of 
the basis states, it is easy to see that if there are more 
ground states than excited states, N g > N e , then the 
dark states will be composed of g-states and some pop- 
ulation will be trapped there. This will remain unaf- 
fected by the radiation; population cannot be removed 
from them using this particular linkage pattern. 

Figure 2 illustrates this accounting procedure. The 
top frame (a) shows a general coupling scheme for the 
sequence J = 2 1 <-> 2. The MS transformation on the 
g — e pump transition produces the description shown in 
the bottom frame (b). In the g set, with this transformed 
basis, there occur two sublevels that have no connection 
with any excited states. Population cannot be removed 
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(a) 




FIG. 2: (Color Online) This sketch shows that when the num- 
ber of ground-state-sublevels exceeds the number of excited- 
state-sublevels then it is impossible to transfer all the popu- 
lation from the ground-state manifold to the final-state man- 
ifold. Frame (a) shows the original coupling configuration. 
Frame (b) shows the couplings after a MS transformation on 
the g — e transition. The empty rectangles represent uncou- 
pled states. This transformation is independent of time, be- 
cause all elements of the coupling share a common time de- 
pendence, p(t). The presence of uncoupled sublevels in the g 
set prevents removal of population from these states; hence 
it is not possible to remove all population from all of the 
ground-state sublevels. 



from these as long as the couplings are those shown in 
the top frame. 

It is easy to see that, had there been more sublevels 
in the e set, such that N g < N e , then every one of the 
transformed states from the g set would be linked to some 
excited state, with consequent possibility for population 
removal. There will also be uncoupled states in the e 
manifold but they are unpopulated and do not affect the 
population transfer. 

The introduction of MS basis states in this way makes 
the g — e linkage pattern quite simple, but by introduc- 
ing a new basis the e — / couplings become more compli- 
cated: generally there will be a connection between each 
transformed e state and each (untransformed) / state, as 
indicated in frame (b). 

Next we consider the e — / coupling. We can repeat 
the previous argument for the g — e coupling with the re- 
placements g <-» e and e <-> f. We obtain, that the Stokes 
field MS transformation yields N < — min{iV e , Nf} inde- 
pendent two level systems for the e — f transition, plus 
\N e — Nf | uncoupled states in the larger one out of the 
e and / sets. It is easy to see that, had there been more 



sublevels in the / set, such that N e < Nf, then every 
one of the transformed states from the e set would be 
linked to some final state, with consequent possibility for 
population removal. Combinig the arguments of the MS 
tarnsormations for the g — e and e — / couplings, we 
obtain that in general, if a non-descending sequence of 
state-degeneracies is fulfilled 



N g < N e < N 



f ■ 



(5) 



then a complete STIRAP-like population transfer from 
the g set to the / set is possible. We emphasize that in 
this case the success of the full transfer is independent of 
the initial state of the system: it can be any pure state 
or a mixed state as well. 

A particularly important special case of degeneracy oc- 
curs when there are dark states but they are insufficient 
to produce complete population transfer. This occurs 
when N g + Nf > N e , but Nf < N e . For example, in the 
linkage ofJ=1^2<-^l there is 1 dark state. Figure 11 
illustrates this situation. 



IV. THE STOKES-FIELD MS 
TRANSFORMATION 



In this section we determine the dark states of the 
Hamiltonian of Eq. (1); these are the adiabatic states 
that will be utilized for the desired adiabatic population 
transfer. In order to simplify the structure of the Hamil- 
tonian, we perform a MS transformation; here we take 
that to be on the e — / couplings (those of the Stokes 
field). In our case the time-independent transformation 
matrix U is defined as 



U 



10 
B 
A 



(6) 



In the top-left corner there is a unit matrix I of dimension 
N g x N g . This leaves the g set of states unaltered. The 
Nf x Nf unitary matrix A transforms the sublevels in 
the final-state manifold. Similarly, the N e x N e unitary 
matrix B transforms the sublevels in the excited-state 
manifold. The constant matrices A and B are defined 
[26] such that by transforming the Hamiltonian Eq. (1) 
with the matrix U through the relation 



UH{t)U^ = 



p(t)P ^ 
p(t)Pt RA s(t)S 
«(t)Sft 



(7) 



we obtain a transformed pump-field coupling matrix 
P = PB\ and a quasi-diagonal Stokes-field coupling 
matrix S = BSA^. By quasi-diagonal we mean that the 
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structure of the matrix is 

'[s 

S x 





(a) 



if Nf > N e , 
if JV> = A e , 

if N f < N e , 



(8) 



where E is a square diagonal matrix with dimension 
N < = min(AT e , Nf). The moduli of the diagonal elements 
are given by the square-roots of the common eigenvalues 
of the Hermitian matrices SS> (of dimension N e x N e ) 
and S' S (of dimension Nf x Nf). The phases of the di- 
agonal elements are obtained by evaluating directly the 
matrix product BSA^ . Some of the diagonal elements 
of S might be zero, meaning that some e — / couplings 
vanish in the MS basis. We here assume that in general 
all diagonal elements of X are non-zero, i.e. it is nonsin- 
gular. We treat in Appendix B the case when this matrix 
is singular. 

In the following subsections we consider the three im- 
portant special cases of degeneracies and derive the adi- 
abatic states of the coupled degenerate systems. 




(b) 




(c) 



A. The case N g < N e < Nf 

We first consider the case when the MS transformation 
on the e — / transition results in Nf — N e > decoupled 
sublevels in the / manifold. The coupling matrix S takes 
the form given in the first row of Eq. (8), and hence the 
Hamiltonian in the MS basis reads 



H(t) = 






p(t)p 








p(*)pt 


HA 


s(i)£ 








s(*)£t 





















e 



FIG. 3: (Color Online) The three stages of the transforma- 
tions, (a) The original coupling scheme, (b) The result of 
the Stokes-field MS transformation, converting the couplings 
between e and / sets into independent one- and two-state sys- 
tems, (c) The result of redefining the states in the g, e, and 
/ sets. 



As with the original RWA Hamiltonian, the only time 
dependence enters through the pulses p(t) and s(t). 

We can treat the system in the same way when Nf = 
N e . Then the coupling matrix S is given by the second 
row of Eq. (8), and we have to omit all zero rows and 
columns from the Hamiltonian of Eq. (9) . In either cases 
the sub-matrix P has dimensions N g x N e , while the 
square matrices X and A have dimensions N e x N e . 

To find the adiabatic eigenvectors 3>fc(t) of H(t) we 
take their elements to have the form 



~x k (t) ' 


9 


Vk{t) 


e 


z k {t) 


f 




r 



where /' denotes the subspace of uncoupled states in the 
/ set. Because these are unlinked to the e set they meet 



the definition of dark states. Their population, if ini- 
tially present, is preserved throughout the time evolu- 
tion. When Nf — N e we simply omit the fourth row 
from this vector (the /' states), i.e. we do not have z' k . 
In Eq. (10) there is no tilde on the x components because, 
unlike the y and z components, these do not transform 
in the Stokes field MS transformation. In Sec. IV B and 
C the x components undergo a MS transformation, as is 
indicated there by a tilde. 

The eigenvectors satisfy the eigenvalue equation 

H(t)# fc (t) = £fc(t)*fc(t) . (11) 

By substituting the Hamiltonian of Eq. (9) and the pa- 
rameterization (10) of the eigenvectors into this equation 
we obtain four sets of coupled linear equations for x k , Vk, 
Zk, and z' k . The solution of these equations provide the 
dark and bright eigenvectors 3>fc(i) defined by Eq. (11). 
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Let us assume that there exists an eigenvalue zero, 
£o = . This is always possible to ensure, by suitable 
choice of the phases of the rotating wave approximation 
and the zero-point of energy. If we can find a solution 
of the eigenvalue-equation (11) for this case, then our 
assumption eo = holds, since the solution of the lin- 
ear equations is unique. After some algebra one can ob- 
tain N g different vectors *&o\t) ,1 — l...N g , that are 
linearly independent of each other, and can make these 
orthonormal 



o 





(0 
o 



(12) 



where N~Q l \t) is a (time dependent) normalization factor. 
Here we have assumed that the matrix S is nonsingular. 
We will discuss separately, in Appendix B, the situation 
when £ is singular. Since the y component of these vec- 
tors is zero, they have no component in the e set; they 
correspond to dark states. To make the dark eigenvectors 
of Eq. (12) orthogonal we require that 

(13) 

for 1 < k < I < N g . The time-dependence of the enve- 
lope functions s(t) and p(t) is arbitrary, and therefore we 
require that the two terms on the left-hand-side (lhs) of 
Eq. (13) be identically zero. The eigenvectors of a Her- 
mitian matrix can be chosen so that they are orthogonal 
to each-other, and therefore the first term on the lhs of 
Eq. (13) is automatically zero. It follows that the vectors 
Xq ^ are the eigenvectors of the Hermitian matrix 



M = P(55' t )" 1 P t = PS-^S-ipt . 



(14) 



There is another set of dark eigenvectors for Nf > N e . 
These follow from the discussion after Eq. (10) and are 
given by 







z /(0 



l = N g + l,...Nf-N e + N g , (15) 



where z' ^ are constant orthonormal unit vectors. These 
dark eigenvectors are clearly orthogonal to those of 
Eq. (12). 

We show in Appendix C that the coupling sequence 
g <-> e <-> / can be rendered to independent three-state 
chains by a suitable hoice of the basis states in the g, e, 
and / sets. Figure 3 illustrates the sequence of trans- 
formations that leads to the construction of the dark- 
state eigenvectors Eq. (12). Frame (a) shows the original 
system, with some couplings. Frame (b) depicts the re- 
sults of the Stokes-field MS transformation of the e and 
/ states. Frame (c) shows the result of the redefinition 



of the g, e, and / sets of states according to Eq. (CI), 
with the resulting set of independent chains. 

The matrix of Eq. (14) may have zero eigenvalues as 
well. If so, the corresponding eigenvectors X ^ satisfy 
the equation 



Pt 



x 



(k) 



0, 



(16) 



since we have assumed that the matrix X is nonsingular. 
Note that here P^ is expressed in the original atomic 
basis. The ith row of the matrix P* describes the cou- 
pling between state i from the e set and the sublevels of 
the g set. The rows of the coupling matrix can be con- 
sidered as vectors that span a subspace of states from 
the g set. The dimension of this subspace is the number 
of linearly independent rows of P\ say Np. Obviously 
we have Np < N e and Np < N g . Therefore, there are 
N g — Np different, nontrivial solutions of Eq. (16). These 
nontrivial solutions provide states that are unaffected by 
the pump field. 

If N g — Np < N e then such an uncoupled state does 

not exist, and the vectors {x^}, k = 1 . . . N g span the 
total g-set manifold. Therefore by choosing a counterin- 
tuitive pulse-sequence for the pump and Stokes pulses, we 
can cause complete transfer of population from the g set 
to the / set by means of independent STIRAP processes. 
For such population transfer to succeed, the conditions 
of the adiabatic evolution should be fulfilled, as we will 
discuss in Sec. V. The success of such population transfer 
is independent of the initial state of the system. It can 
be any single state, an arbitrary coherent superposition 
of states or even a mixed state, see Sec. V. 

If Np < N g then some g-set sublevels are decoupled 
from the pump field, hence in general it is then impossible 
to move all the population from the g set. Part of it is 
trapped in dark states. 

The other 2N e adiabatic eigenvectors belong to non- 
zero eigenvalues. They can be obtained in the form 



**(*) - 



A4(i) 



p{t)Pm{t) 

£k{t)Vk{t) 
«(t)Sty fc (t) 




k = 1...2N e 



where A4 (t) is a normalization factor and y k (t) satisfies 
the eigenvalue equation 



p(i) 2 ptp + w (t) 2 ££tJ yk( fi = £k (t)[e k (t) - hA]y k (t) . 

(18) 

Because they contain component states from the e set, 
these are bright states. Although for population transfer 
we use the dark states of Eq. (12), we need the bright 
states to find the adiabaticity conditions; see Sec. V. 

In summary: in this subsection we have shown that 
when N g < N e < Nf, under very general conditions the 
complete population from the g set can be transferred 
to the / set of states. Once we have fixed the pulse- 
shapes, polarizations and phases, complete transfer can 
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be obtained for any arbitrary initial state from the g set. 
The eigenvectors k = 0. . . 2N e in the original bare 

atomic basis can be obtained as 



*fc(t) 



Af k (t) 



" x k (t) ' 




Vk(t) 


1 


Zk(t) 


N k (t) 








x k (t) 



At 



Zk(t) 





(19) 

Moreover, with this method it is possible not only to 
transfer populations, but to create superposition states 
in the / set. We will consider this possibility in Sec. VI. 



B. The case N g > N e > N f 

According to the considerations presented in the begin- 
ning of Sec. Ill, we cannot expect that all the population 
from the g set can be removed when N g > N e > Nf. 
However, a part of the population can be removed and 
with this we can create coherent superposition states in 
the / set. In order to find the dark- and bright states of 
the system we proceed in the same way as in Sec. IV A, 
but now with the MS transformation involving the pump 
transition 



U 



BOO 
A 
1 



(20) 



We look for the eigenvectors of the transformed Hamil- 
tonian in the form 



**(*) - 



" x k (t) ' 


9 


S'fc(*) 


9' 


Vk(t) 


e 


. z k {t) _ 


/ 



(21) 



The vectors x' k (t) describe the population in those states 
of the g set that are decoupled from the pump field. 
There are N g — N e dark states in the g manifold, and 
these can be written in the form 



*?>(*) 





-'(0 
"o 






(*) 



(22) 



where the vectors {x , (i)} form an orthonormal set. 
The population cannot be removed from these states. 
The rest of the dark states are obtained in the manner 
used for Eq. (12). They can be written as 








(fe) 




(23) 



= BPA\ with II a diagonal coupling ma- 
trix of dimension N e x N e , and S = AS. We require 



where 



n 

o 



orthogonality for the dark states Eq. (23). Hence the 
constant vectors Zq are chosen so that they are eigen- 
states of the Hermitian matrix S , tn _1 nt _1 S, in direct 

(k) 

analogy with the way the constant vectors x y were cho- 
sen earlier in Sec. IV A. 



C. The case N g ,N f < N e 



Here we consider the situation N g ,Nf < N e . We 
will show that under these conditions the dark states of 
the system can be identified by means of two sequen- 
tial MS transformations. The first MS transformation is 
performed among the e and / sets of the Stokes transi- 
tion, as in subsection IV A. The transformation matrix 
is given by Eq. (6). As a result, the coupling matrix S 
of the Hamiltonian (1) takes the quasi-diagonal form of 
the third row of Eq. (8) . Therefore, the Hamiltonian in 
the MS basis reads 



H(t) = 



p{t)P a p(t)P b 
p(t)Pt HA s(t)E 
p{t)Pl ] HA 

s(t)£t o 



(24) 



The diagonal square matrix S has dimension Nf x Nf. It 
can be readily seen that there are N e — Nf states in the e 
set that are not coupled to the / set. We call these uncou- 
pled levels, whereas the other subset of coupled excited- 
state-sublevels are called active. In the Hamiltonian of 
Eq. (24) the pump coupling matrix is partitioned into 
two sub-matrices: the matrix P a of dimension N g x Nf 
describes couplings between the g set and the active MS 
states of the e set. The other sub-matrix of dimension 
N g x (N e — Nf) is associated with the transitions between 
the states of the g set and the uncoupled states of the e 
set. The result of this transformation is illustrated in 
Fig. (4b). As the figure shows, we cannot identify clearly 
the dark states, because in general all the states of the g 
set are coupled to all of the e set. Therefore, we perform 
a second MS transformation, involving the g set and just 
those states of the e set that are decoupled from the / 
set - two in the present example. The result is illustrated 
in Fig. (4c). In this example there is one (/-set state that 
couples solely to an active MS state of the e set because 
the other two have, by means of the MS transformation, 
been linked to the two uncoupled e states. The popula- 
tion can be moved from this g state to an / state. The 
middle e state is coupled to all three g states. Conse- 
quently, if the two spectator g states are populated, they 
disturb the complete population transfer from the middle 
g state into / states, and the population transfer process 
will place population into the e and / states. 

In general, the transformation matrix of the second MS 
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(a) 




transforms the uncoupled states of the e set. The two 
unit matrices J are of dimension Nf x Nf. For N e — Nf < 
N g the transformation yields 



A'P a = 

A'P b B n = 



P 
P' 


n 



(26a) 
(26b) 



(b) 



(c) 




FIG. 4: (Color Online) The three frames show the stages 
of the MS transformations when the sequence of degenerate 
states violates the condition Eq. (5): Frame (a) depicts the 
original coupling scheme , here Nf < N e . Frame (b) shows 
the result of the Stokes field MS transformation, converting 
the couplings between e and / sets into independent one- and 
two-state systems. The one-state systems are in the e set. 
Frame (c) shows the result of the pump field MS transforma- 
tion followed by the redefinition of the sublevels in the g, e, 
and / sets according to Eq. (CI), leading to one dark state. 
This dark state is associated with the middle three-state link- 
age, indicated by heavy lines. In addition, the middle e state 
is coupled not only to a single g state but to the two others 
as well. As a result, the populations in the two spectator g 
states may disturb the complete population transfer from the 
middle g state, see text. 

transformation is defined as 

" A' 0" 

, J ( 

U B' ■ [l *> 

.0 /. 

The N g x N g unitary matrix A' transforms the g set, 
whereas the (N e — Nf) x (N e — Nf) unitary matrix B' 



where the matrix P is of dimension [N g — (N e — Nf)} xNf, 
P' is of dimension (N e — Nf) x Nf, and II is a diagonal 
matrix of dimension (N e -N f )x(N e -N f ). For N e -N f = 
N g we find 

A'P a = P', (27a) 
A'P b B n = fl, (27b) 

where the matrix P' is of dimension N g x Nf and II is 
of dimension N g x N g . We do not have P in this case. 
Finally, for N e - Nf > N g we get 

A'P a = P' , (28a) 
A'P b B'^ = [no], (28b) 

where the matrix P' is of dimension N g x Nf and II 
is of dimension N g x N g . Just as with the conditions 
N e — Nf — Ng, we do not have a matrix P in the present 
case either. 

In general, none of the diagonal elements of the matrix 
II are zero. Therefore, when N e — Nf > N g there are 
no states in the g set that are coupled solely to active 
MS states in the e set. It follows that no dark state 
can be identified in the system and hence a STIRAP-like 
population transfer is impossible. 

In special (but important) cases it may occur that some 
diagonal elements of II vanish. Then the system has 
such MS g-set states that are coupled only to active MS 
states of the e set, hence the system has dark states and 
a STIRAP process is possible. We will reconsider this 
case later in this subsection. 

In all three cases N e - Nf > N g , N e ~ Nf = N g , and 
N e — Nf < N g some initial population of the g set cannot 
be included in the dark states of the system in the general 
case of arbitrary initial superposition of g states. We 
conclude that in general in the case of N g , Nf < N e it is 
impossible to remove all the population from the g set in 
a STIRAP-like population transfer process. Exceptions 
occur when the matrix II is identically zero. Then the 
uncoupled MS states of the e set are decoupled not only 
from the / set, but also from the g set. 

When N e — Nf < N g the second MS transformation 



9 



produces from the Hamiltonian of Eq. (24) the matrix 



H'(t) = 



p(t)P 
p{t)P' p(t)U 
hA 
hA 
s(t)£t o 






s(t)E 








p{t)pt p(t)p' f 
o P (t)W 



" (29) 

The situation N e — Nf > N g can be treated similarly. 
In order to find the dark states of the Hamiltonian of 
Eq. (29) we proceed in the same way as in Sec. IV A. 
The eigenvectors are parameterized as 



(30) 



The eigenvalue equation is defined by Eq. (11). By in- 
serting the Hamiltonian of (29) and the parameterization 
of the eigenvector Eq. (30) into Eq. (11) we obtain five 
groups of coupled linear equations for x k , x' k , y k , y' k , 
and z k . The No linearly independent dark eigenvectors 
3>Q ^ (t) that belong to the eigenvalue e 





" x k (t) ' 


9 




3b(*) 


9' 


**(*) = 


Vk{t) 


e 






e! 




. z k (t) . 


f 



0, are given by 



JVo(t) 



s(t)x 
s(t)x' 







.p4; 



(31) 



where 
condition 



the constant vector x' should satisfy the extra 



Ila^ = . 



(32) 



This condition says that in a dark state no population can 
be in those g-set states that are linked to uncoupled e-set 
states. An example to this configuration is shown later 
in Sec. VIC. The dimension No of the dark subspace is 
equal to N g — (N e — Nf) plus the dimension of the zero- 
subspace of the matrix II, Eq. (32), where we assumed 
that N g > (N e -Nf). For N g < (N e -N f ) the dimension 
of the dark subspace is equal to the dimension of the zero- 
subspace of the matrix II. 

It is useful to orthogonalize the dark states of Eq. (31). 
The orthogonality relation is given by Eq. (13). In this 



(0 ™'(0lT 



= l...Nr 



case we find that the vectors [xq x q 
should be the eigenvectors of the Hermitian matrix 



P 

P' 



S -it S -i[ptp't] ? 



(33) 



with the restriction of Eq. (32). 

The other eigenvectors, belonging to non-zero eigen- 
values, are given by 



*fc(t) = 



A4(i) 



p(t)Py k {t^ 
P(l) P'yiAh ■ Ilyti/i 

e k (t)y k (t) " - ( 34 ) 

£fc(*)Vfc(*) 
a(t)Sty fc ( t ) 



where Afk(t) is a normalization factor and the vector 
[yk(t) y'k(t)] T satisfies the eigen-equation 



P (t) 2 (ptp + p ,] p'1+ s(i) 2 sst P (t) 2 p'tn 
P {t) 2 WP' p(t) 2 itfn 



Vk(t) 



= Sk{t)[s k (t) - hA] 



Vk{t) 

m 



(35) 



The states of Eq. (34) are bright states, because they 
include components from the e set. 

The eigenvectors of the Hamiltonian H(t) in the bare 
atomic basis can be obtained as 



x k (t) 
Vk{t) 

W) 

z k (t) 



Pt 



x k (t) 
x' k (t) 

Vk(t) 
B'^y'S) 



A^z k (t) 



(36) 



V. ADIABATICITY CONDITIONS AND TIME 
EVOLUTION 

In Sec. Ill we have presented the dark states of our 
degenerate system. Once we have the dark states we 
may consider adiabatic evolution of the system in the 
dark subspace. There are two questions that should be 
addressed in connection with adiabatic evolution: 

(1) What are the conditions needed to ensure adiabatic 
evolution? 

(2) If there are several degenerate dark states of a sys- 
tem, in general there are nonadiabatic couplings 
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among them. How can we find the time evolution 
of the system in this case? 

To answer the first question we apply the basic theory 
of adiabatic evolution [29], which assures that the evo- 
lution is adiabatic if any nonadiabatic couplings among 
the adiabatic states are negligible compared with their 
energy separation. In our model system we have a dark 
subspace that is spanned by states that have eigenvalue 
zero. The other adiabatic states, the bright states, have 
non-zero eigenvalues. Because we want the state vector 
to remain in the dark subspace, we require that the dark 
subspace be separated from the bright one, as expressed 
by the condition 



t»\(*W(t)\* k (t))\<Z\e k (t)\, 



(37) 



where I = 1 . . . No and k = 1 . . . Nb, with Nb being the 
number of bright states. The dot denotes time derivative. 
We may insert into Eq. (37) any set of dark and bright 
states from Sec. III. For example using the dark states 
given by Eq. (12) and the bright states from Eq. (17) we 
find 



s{t)p{t) - p(t)s{t)l (x^\P\y k (t)) «|e fc (t)| 

(38) 

This formula closely resembles the adiabaticity condition 
of the conventional nondegenerate three-level STIRAP 

[3], 



S(t)P(t) - P(t)S(t) 



sin if 

COS(/5 



cot If 
tany 



with flo(t) and ip defined as 



fto(t) = y/P 2 (t) + S 2 (t), tan2^ = 



fio(t) 



(39) 



(40) 



respectively. The second factor on the lhs of Eq. (38) 
contains the time-derivatives of the envelope functions of 
the pump and Stokes pulses, whereas the third factor in- 
volves the element of the matrix P between the Ith dark 
state at the initial time and the excited-state-amplitudes 
of the fcth bright state at time t. On the rhs e k is the 
eigenenergy associated with the fcth bright state. When- 
ever the adiabaticity conditions Eq. (38) are fulfilled for 
all dark and bright states of the system, then the dark 
and bright subspaces evolve independently. 

There remains the task of determining the time evo- 
lution in the dark subspace. When there are several de- 
generate dark states there are usually nonadiabatic cou- 
plings among them. In case of the tripod system [17, 18], 
due to the special choice for the time-dependence of the 
Stokes pulses, the two dark states mix throughout the 
population transfer process. If the dimension of the dark 
subspace is larger than two, then in general there is no ex- 
act analytic solution [21, 22]. In our case the situation is 
much simpler. The nonadiabatic coupling between a pair 



of dark states is ($q (i)|$ (t)} . By evaluating this ex- 
pression for any pair of dark states from Sec. IV we always 
get identically zero. Hence the dark states do not mix 
throughout the whole transfer process. This property 
simplifies the calculations considerably, since the time 
evolution operator in the dark subspace is given by 



N D 



(41) 



i=i 



For example, for N g < N e < Nf the vectors <J>Q^(t ) are 
equal to a;^ from Eqs. (12) and (13). Once we define the 
density matrix g(to) of the system at time to the density 
matrix at any later time is given by 



Q{t)=U{t,to)g(to)U^(t,t ) 



(42) 



provided that the adiabaticity conditions Eq. (38) are 
satisfied. Note that this formula is valid only if f?(io) lies 
entirely in the dark subspace. It follows from Eqs. (41) 
and (42) that any pure or mixed initial state of the system 
occupying the dark subspace of the g set is transferred to 
the / set in the course of the population transfer process. 
In the case of a pure initial state ^(to) we have 



*(i) = l7(t,to)*(to). 



(43) 



Here again, \&(io) must have components solely in the 
dark subspace. (Were the state vector to have compo- 
nents initially in the bright subspace, then adiabatic evo- 
lution would maintain such presence. Because the excited 
states undergo spontaneous emission, their populations 
have the potential to interrupt the coherence of the dy- 
namics and thereby to diminish the population transfer.) 



VI. SOME EXAMPLES 

In this section we demonstrate through some exam- 
ples the usage of our method. To be specific, we consider 
atomic transitions where the origin of the degeneracy is 
the set of degenerate magnetic sublevels of angular mo- 
mentum states in the absence of a magnetic field. Our 
purpose is to present some typical configurations that 
may occur in realistic situations. 



A. The J = 1 <-> 2 <-> 3 linkage 

In this example we consider the linkage J = 1 «-> 2 <-» 
3, shown in Fig. 5 and assume that only fields are 
present. In this case there are two independent coupled 
systems: the one with M g = 0, M e = ±1, and Mf 



0,±2 (shown as dashed lines); and the other one with 
M g = ±1, M e = 0, ±2, and M f = ±1, ±2 (shown as full 
lines). The first one has been studied in ref. [25], hence 
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(a) 




(b) 




(c) 



t r 



1 / 
1 / 
i / 
i / 

.V. 



t r 



/ 



8 



FIG. 5: (Color Online) The coupling configuration 
for the J = 1 <-> 2 <-» 3 linkage with only a polarized 
coupling fields. The system separates into two independent 
subsystems; the smaller one is shown with dashed lines, the 
larger one with solid lines. Frame (b) shows the result of the 

Stokes-field MS transformation. Frame (c) shows the 
redefinition of the states in the g, e, and / sets according to 
Eq. (CI). 



we do not consider it here. For the second, larger system, 
the pump coupling matrix P is given by 



h l 
27! 



p Ve p 
TE^ 





(+) 



(44) 



whereas the Stokes coupling matrix reads 



(+) 



(-) 



n 






(+) 



(45) 




4 +) 


= !lpe* 


cos 77 , 


n<r> 


= !lpe* 


sin 77 , 




= ft s e^ s 


cos 8 , 


4-) 


= risers 


sin 8 , 



FIG. 6: (Color Online) The eigenvalues of the matrix SS^ , 
Eq. (45), as a function of the polarization of the Stokes field. 
The eigenvalues are measured in the units of (HQs) 2 ■ 



The numeric factors in front of the fl-s describe the 
Clebsch-Gordan coefficients. The Rabi frequencies are 
parameterized as 



(46a) 
(46b) 
(46c) 
(46d) 



where the amplitudes f2p,s are nonnegative. The angles 
rj and 8 characterize the pump and Stokes field polariza- 
tions, respectively. 

Here N g < N e < Nf (2 < 3 < 4), hence the derivation 
in Sec IV A can be applied. As a first step, we have to 
perform the Stokes field MS transformation. The eigen- 
values Afc of the matrix SS^ are given by the roots of a 
cubic equation, see Eq. (Dl). We display them in Fig. 6 
as a function of the polarization angle 8. They are never 
zero, hence the complete adiabatic population transfer is 
possible for any polarization of the Stokes field. How- 
ever, their amplitudes depend on the polarization, which 
affects the adiabaticity conditions Eqs. (37) and (38). 
The Stokes field MS transformation matrices A and B, 
Eq. (6), can be calculated in a straightforward manner; 
they are shown in the Appendix D. Since Nf = N e + 1, 
the Stokes field MS transformation yields a transformed 
coupling matrix S in the form of the first row in Eq. (8). 
The diagonal part S is given by 



20 



A, 









/A2 





(47) 



There are 2 + 4 — 3 = 3 dark states in this system: one 
is in the / set, an uncoupled state. The space of g is 
two-dimensional, N g — 2, and hence there are two dark- 
states in the form of Eq. (12). The vectors associ- 
ated with these two dark states are the eigenvectors of 
the Hermitian matrix M of Eq. (14) and are given in the 
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FIG. 7: (Color Online) Upper frame: the pulse sequence used 
for the population transfer process in the coupled angular 
momentum system J = 1 <-> 2 <-> 3. Initially only the state 
\g, J g — l,M g = 1} was populated. Lower frame: The popu- 
lation evolution. 



Appendix D. The two dark states are obtained by insert- 
ing their structure into Eq. (12) or, for the bare atomic 
basis, into Eq. (19). 

We have performed numerical simulations to check the 
validity of our analytic results. In Fig. 7 initially the sys- 
tem was in the state \g, J g = 1, M g = 1). The envelope 
functions of the pump and Stokes pulses, respectively, are 
P (t) = exp(-[t - 3] 2 /6 2 ) and s(i) = exp(-[t + 3] 2 /6 2 ); 
and flp = 52, f2g = 42. The polarizations are character- 
ized by 77 = 1.3376 rad, and 8 = 0.4636 rad. The phases 
are chosen randomly as <pp = 1.1814 rad, ipp = rad, 
<f) S = 1.8925 rad, and ip s = 2.8198 rad. The detuning 
A is set to zero. We have found again very good agree- 
ment between the analytic calculations and the numeric 
simulation. 

We also considered a mixed initial state, when the ini- 
tial state of the system is chosen as half of the population 
is placed on each of the \g, J g — 1, M g = ±1) states, and 
the coherence is zero between them. The numerically 
calculated dynamics is shown in Fig. 8. We can see that 
despite of the mixed initial state, the complete popula- 
tion can be transferred from the g set to the / set. The 
pulse sequence is the same as in the previous example. 







\p(g-l) 






P(g,l)\ 












P(f,3) 








Pi f,- 1 i 








P(f,l) 








P(f,-3) 



t farb. unitsl 

FIG. 8: (Color Online) Same as Fig. 7, but for a mixed initial 
state, Pi(l, —1) = P:(l, 1) = 1/2, all initial coherences are 
zero. 



B. The J = 1 



1 linkage 



As another example we consider the linkage J = 1 
1 <-> 1 shown in Fig. 9. In this case N g = N e = Nf, and 
hence the derivation in Sec. IV A is applicable. This is a 
counter-example to the general condition of Eq. (5): even 
though the condition Eq. (5) is satisfied, in this case the 
complete removal of an arbitrary population distribution 
from the g set is impossible in the STIRAP way. 

The coupling matrices S and P in the Hamiltonian 
Eq. (1) are given by 



X 



2^6 



(t) 



-n x 



-n 



(+) 

X 







(-) 

X 



-n 



(+» 







(-) 

X 



X 

(tt) 
X 



(48) 



for X = S or P. The factor l/y/6 and the ± signs 
describe the Clebsch-Gordan coefficients. The Rabi fre- 
quencies f2^' correspond to the er + , a~ , and tt polar- 
izations, respectively. Note that a selection rule nullifies 
transitions M = <-> M = 0. 

As described in Sec. IV, we perform the Stokes-field 
MS transformation to diagonalize the Stokes coupling 
matrix S. The eigenvalues of the matrix SS' provide 
the squared moduli of the diagonal elements of the ma- 
trix S = S = BSA\ Eq. (7). They are given by 



0, ±- m 



10 



,(rms) 



(49) 



with ( £1 



(rms] 



in 



(+)|2 



(")|2 



fi^l 2 . One of 



the eigenvalues is always zero and therefore, although 
the system satisfies the condition for complete popula- 
tion transfer, Eq. (5), the null Rabi frequency prevents 
complete transfer. 

Fig. 10 demonstrates the population transfer in this 
system. Initially the system was in the state (\g,J g — 



l,Mn 



\g,Jg = l,Mg = 0) + \g,Jg = l,Mg 
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(a) 




(b) 




(c) 
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P(f,0) 
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FIG. 10: (Color Online) Upper frame: the pulse sequence 
used for the population transfer process in the coupled angu- 
lar momentum system J = 1 <-» 1 <-> 1. Initially the state 
(| ff) J g = 1, M g = -1) - Iff, J g = 1, M 9 = 0} + |ff, J 9 = 1, M g = 
l})/%/3 was populated. Lower frame: The population evolu- 
tion. Part of the population is left in the g set, because some 
of the MS Rabi frequencies vanish. 



FIG. 9: (Color Online) Same as Fig. 5 for equal state- 
degeneracies. For equal J-s, the M = <-+ transition is 
dipole-forbidden, and we cannot select a basis such that there 
occur couplings between all pairs of states of the degenerate 
sets. This restriction results from the property of the Clebsch- 
Gordan coefficients. Frame (c) shows the redefinition of the 
states in the g, e, and / sets according to Eq. (CI), leading 
to two independent three-state linkages. 



C. The J = 1 <-» 2 <-> 1 linkage 

In our last example we consider the linkage J = 1 
2 <-» 1, shown in Fig. 11 and assume that only a fields 
are present. In this case there are two independent cou- 
pled systems: the one with M g — 0, M e = ±1, and 
Mf = (shown as dashed lines) discussed recently in 
ref. [30], and the other one with M g = ±1, M e = 0, ±2, 
and Mf = ±1 (shown as full lines). This is a twin di- 
amond configuration. For the larger system, the pump 
coupling matrix P is given by 



l))/v3. The envelope functions of the pump and Stokes 
pulses, respectively, are chosen as pit) = exp(— [t — 
2] 2 /4 2 ) and s(t) = exp(-[i + 2] 2 /4 2 ), and Q P = fl s = 30. 
The intensity is equally distributed among the <j + , er~, 
and 7r components of the exciting fields. The detuning 
A is set to zero. We have found excellent agreement 
between the analytic calculations and the numeric sim- 
ulation. The adiabaticity conditions, Eq. (38), are also 
fulfilled throughout the relevant part of the population 
transfer process. 



73 p 







(50) 



whereas the Stokes coupling matrix reads 



f1~+, 





5 "S 



( + ) 

S 

(") 



(51) 
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The parameterization of the Rabi frequencies * s 
given by Eq. (46). Here N g ,N f < N e (2,2 < 3), hence 
the derivation in Sec IV C is applicable. The sequence 
of the dimension of the subspaces violate the condition 
N g < N e < Nf, therefore, in general a STIRAP-like 
complete population transfer is not possible. However, 
this is another counter-example to the general condition 
of Eq. (5): even though the condition Eq. (5) is vio- 
lated, we show that the complete removal of an arbitrary 
population distribution from the g set is possible in the 
STIRAP way for a special choice of pulse polarizations 
and phases. 

As usual, we start with the Stokes field MS transfor- 
mation. The eigenvalues of the matrix SS^ are given by 
the roots of a quadratic equation, which read 

7 1 



Al ' 2= 100 ± iooV24cos 2 2(9+l. (52) 
The Stokes field MS transformation matrices A and B, 



Eq. (6), can be calculated in a straight forward manner, 
they are shown in the Appendix E. Since Nf = N e — 1, 
the Stokes field MS transformation yields a transformed 
coupling matrix S in the form of the last row in Eq. (8). 
The diagonal part £ is given by 



AT 
o 



A2 



(53) 



The eigenvalues of Eq. (52) are always positive, hence 
this matrix is nonsingular for any polarization of the 
Stokes field. The Stokes field MS transformation yields 
two e — / linkages and an e state which is not coupled to 
any / state, see Fig. lib. Now, following the derivation 
of Sec IV C we perform a second MS transformation for 
the pump field. The transformation matrix is given by 
Eq. (25). In our case, the 2x2 unitary matrix A' is 
defined as 



A' 



cosf 



e -*«>s-0s) sin( 



■hi(<t>s-^s+4>p +4>p ) sin0 — e~'5 i ('/'s-0s+0-p+i/'-p) 



COS) 



(54) 



while the matrix B' is a scalar now, and chosen as unity. Since in this case N e — Nf < N g (3 — 2 < 2), the transformed 
pump field coupling matrix takes the form of Eq. (26). The matrix II is a scalar, that reads 



n 



2V3V2 — cos 2 29 



(cosr/cos6»elH'/>s-0s+<^- l />p) - sin^sin^-^ 3 -^"^-^)) . 



(55) 



This is nonzero in general, hence one of the g states is 
linked to the uncoupled e state. Therefore, there is one 
dark state in the system, which reads 



N^(t) 



sit) 







• (56) 



This dark state is associated with the three-state linkage 
in the middle of Fig. 11c, indicated by heavy lines. 
However, for 



i>s - (t>s + 4>p - i>p = k-K , (57a) 
0+{-l) k V = l*, (57b) 



where k is an integer, the scalar II vanishes. As a result, 
the uncoupled e state becomes decoupled from the g state 
as well. Therefore, beside the dark state of Eq. (56) there 



is an other one 



*?>(*) = 





s(t) 




-p{t)ll- l B a P^A^ 



(58) 



In summary: complete population transfer is possible 
from the g set to the / set for the special choice of pulse 
polarizations and phases Eq. (57). It is important to note 
that the condition for complete transfer Eq. (57) is equiv- 
alent to that for the diamond configuration [30]. Hence, 
the population from the total g set can be transferred 
into the / set if the condition Eq. (57) is fulfilled. 

Fig. 12 demonstrates the population transfer in the 
twin diamond configuration. Initially the system was in 



the state cos(a)\g, J g — l,M g 



-1) + sin(a)|g, J g = 



1, M g = 1) with a — arctan(l/3). The envelope functions 
of the pump and Stokes pulses are chosen as in Sec. VI B. 
The polarization of the pump and Stokes pulses were cho- 
sen as n — 2ir/5 and 9 = —tt/7, respectively. All phases 
of the pulses are zero. The detuning A is set to zero. 
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(a) 



(b) 



(c) 




FIG. 11: (Color Online) Same as Fig. 5 for the J = 1 <-> 2 <-> 1 
linkage with only a ± polarized coupling fields. The system 
separates into two independent subsystems; The smaller one 
is shown with dashed lines, the larger one with solid lines. 
Frame (b) shows the result of the Stokes-field MS transfor- 
mation. Frame (c) shows the result of the pump field MS 
transformation and redefinition of the states in the g, e, and 
/ sets according to Eq. (CI). In general, there is one dark 
state in the larger system which is associated with the three- 
state linkage indicated by heavy lines. 

After the pulse sequence has passed, some population is 
left in the g and e sets because the polarizations of the 
pulses violate the special condition for complete trans- 
fer Eq. (57). Finally, in Fig. 13 the polarizations of the 
pump and Stokes pulses are chosen so that the special 
condition Eq. (57) is fulfilled. Then, a complete popu- 
lation transfer occurs, all population from the g set is 
moved into the / set. 



VII. SUMMARY 

We have considered the extension of the well-known 
STIRAP process in degenerate systems in which N g de- 
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FIG. 12: (Color Online) Upper frame: the pulse sequence 
used for the population transfer process in the coupled angu- 
lar momentum system J = 1 <-> 2 <-> 1. Initially the state 
cos(a)\g, J g = 1, M g = — 1) +sm(a)\g,Jg = l,M g = 1} with 
a = arctan(l/3) was populated. Lower frame: The popu- 
lation evolution. The pulses are chosen so that the special 
condition for complete transfer Eq. (57) is violated, hence 
part of the population is left in the g and e sets. 



generate states of the g set are coupled by means of a 
pump pulse to N e degenerate states of the e set, which 
in turn are linked by the Stokes pulse to Nf degenerate 
states of the / set. We have shown that such a generalized 
STIRAP process is always possible if the succession of 
state-degeneracies is nondecreasing, i.e. N g < N e < Nf, 
and the number of non-vanishing MS Rabi frequencies 
is at least N g for both the pump and Stokes couplings. 
When such conditions hold, then for arbitrary couplings 
among states (e.g. arbitrary elliptical polarization of 
electric dipole radiation between magnetic sublevels) it 
is possible to obtain complete adiabatic passage of all 
population from the states of the g set into some combi- 
nation of states of the / set. In this process the initial 
state is arbitrary, it can be any pure or mixed state that 
occupy the g set. 

An important exception from the above rule occurs in 
coupled angular momentum systems, when J g — J e — 
Jf. Then, due to the symmetry of the Clebsch-Gordan 
coefficients some couplings vanish, which results in in- 
complete transfer. 

We have examined the possibility of adiabatic passage 
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FIG. 13: (Color Online) Same as Fig. 12, but now the special 
condition for the pulse polarizations Eq. (57) is fulfilled. Up- 
per frame: the pulse sequence used for the population transfer 
process. The polarizations for the pump and Stokes pulses are 
■q — 2n/5 and 6 = tt/10, respectively. Lower frame: The pop- 
ulation evolution. All population is transferred from the g set 
into the / set. 



when this restriction on degeneracies does not hold. We 
have shown that part of the population can be trans- 
ferred to the / set. We have also pointed out that, for 
certain choices of the polarizations of the coupling fields, 
complete adiabatic population transfer can be obtained. 

We have demonstrated that our scheme can be a pow- 
erful tool for coherent control of the quantum state in 
a degenerate system: in our proposal the selective ad- 
dressing of individual states in the degenerate sets is not 
required. Nevertheless, the final state can be tailored by 
varying the polarizations and the relative phases of the 
coupling fields. We have shown through some specific 
examples that the control of the final superposition state 
is possible; the level of control depends on the system 
under consideration. 
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APPENDIX A: DIPOLE TRANSITION 
MOMENTS 

A common situation where degeneracy occurs is when 
the atomic states are eigenstates of angular momentum, 
bearing the labels J and M. Then the dipole moments 
can be expressed in terms of Clebsch-Gordan coefficients 
and reduced matrix elements. For the pump transition 
(<7 — e) the general pattern of the dipole-transition matrix 
elements, for arbitrary polarization, is 



„ V ^"Jg + 1 1 J 



/'-.;='.'//"! 2..'-7" =rr=r ' 1 ; = ■!...#' ' 

(Al) 

where (g|/Li|e) is the reduced matrix element and pa- 
rameterizes the contribution of spherical component q to 
the interaction. The Stokes transition moments are sim- 
ilarly written as 



lJ) 2^ € i J277TT 1 



i = 1 ■ ■ ■ N e 
J 1 • • • \> 

For vibrational transitions in molecules the reduced 
matrix element must include a Franck- Condon factor. 



APPENDIX B: SINGULAR COUPLING MATRIX 

S 

Let us consider the Hamiltonian Eq. (1) in the MS 
basis, Eq. (7). The MS transformation of the coupling 
matrix S may result in three different forms, shown in 
Eq. (8). We obtain a diagonal matrix X to which are ap- 
pended either rows (if Nf < N e ) or columns (if Nf > N e ) 
of zero values. In the discussions of Sec. Ill we have as- 
sumed that the matrix S is nonsingular. Here we con- 
sider the case when some diagonal elements of X are zero. 
Let us choose the MS transformation matrices A and B 
in Eq. (6) in such a way that the zero diagonal elements 
appear in the bottom right corner of S. This non-zero 
part is denoted by Sc. Let the dimension of this matrix 
be Nq x Nc- In this notation, instead of Eq. (8) we have 



S = 







(Bl) 
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where the number of all zero rows is N e — Nc and the 
number of all zero columns is Nf — Nc- From this form of 
the Stokes coupling matrix it is clearly seen that we have 
N e — Nc uncoupled MS states in the e set and Nf — Nc 
uncoupled MS states in the / set. By inserting the cou- 
pling matrix Eq. (Bl) into the transformed Hamiltonian 
of Eq. (7) and performing a second MS transformation 
as in Sec. IV C among the g set and the uncoupled MS 
states of the e set we get 



H(t) 






p(t)pt p(t)p* 

















p(t)U 


hA 









s(t)£ c 






p(t)p 
p(t)p' 

hA 



s(t)% 


(B2) 

This Hamiltonian is almost identical with the one in 
Eq. (29). The difference is that here on the bottom of the 
matrix we have some rows of zero values as well as some 
columns of zero values to the far right. The adiabatic 
states of this Hamiltonian can be found as in Sec. IV C. 
The eigenvectors are parameterized as 



Vk 

Zk 



(B3) 



The eigenvalue equation yields the set of equations as in 
Sec. IV C plus one more equation for z' k 



= e k z k . 



(B4) 



When looking for the eigenstates belonging to the eigen- 
value zero we set eo = in Eq. (B4). Since z' does not 
appear in the other equations, its value is determined 
from the initial condition of the system. Our usual as- 
sumption is that initially only the states of the g set are 
occupied, therefore, z' = 0. For the eigenstates with 
non-zero eigenvalues the only way to satisfy Eq. (B4) is 
to set z' k to a null vector , z' k = 0. The eigenstates associ- 
ated with non-zero eigenvalues e k are given in Sec. IV C. 



APPENDIX C: LINEARIZATION OF THE 
COUPLINGS g^e^ f 



The construction of the dark state Eq. (12) can be 
understood as follows. We introduce three sets of states, 



defined in the g, e, and / sets, respectively 
g set: tp g 
e set: ip e 



x$' , l = l---N g , 



.(') 
o 

1 pt-C) 



l = l---N„ 



(Cla) 
(Clb) 



+ N e — N g other linearly independent states 



f set: V/ = 



ipU(') 



' 



= l---N n 



(Clc) 



+ Nf — N g other linearly independent states . 

The vectors are orthonormal by construction; 

and JVj 1 ^ are appropriate normalization factors for the 
other components. The states in the g and / sets of 
Eq. (CI) are orthonormal, but the states in the e set of 
Eq. (Clb), though linearly independent and providing a 
complete set of excited states, are not orthogonal. The 
dual counterpart [28] of the e set of Eq. (Clb) reads 



dual e set: ^ e 



^4 )T P^-\ (C2) 
l = l---N„, 



+ N e — N g other linearly 
independent states. 

The vectors of these two sets are mutually orthogonal 



^(0 ~« 



3fcZ • 



(C3) 



In the basis defined by Eqs. (CI) the Hamiltonian of 
Eq. (9) reads 



H(t) = 



p(t)Qi 

p(t)Q 2 ftA «(t)Si 

s(t)S 2 





(C4) 



where Q2 and Si are diagonal matrices with elements 

l = l---N g , (C5a) 



(Q 2 )ii = 

and 
(Si)h = 



Nf 2 



A4 W ' 
1 



(Q 2 )n, l = l---N g , (C5b) 



respectively. It can be verified that the matrix elements 
of pt is zero between the rest of the dual e states and 
the g states 

(il> e |P t |^)=0, k = N g + l---N e , l = l---N g . 

(C6) 

Similarly, the matrix elements of £ is zero between the 
rest of the dual e states and the first N g f states 

(V'e |S|V / ) = 0, k = N g + l---N e , l = l---N g . 

(C7) 
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The matrix elements of the other two symmetric, non- 
diagonal matrices Q\ and £ 2 are given by 

(Qi)ik = j^{4 )T \P&\ X W), (C8a) 
and 

(S 2 ) ife = j^(Qi)ik- (C8b) 

The dark states of the Hamiltonian (C4) can be obtained 
in the same manner as in the above derivation that led 
to the dark states Eq. (12). The population transfer is 



described by the equation 

pWtf^ + a(tjzo =0, (C9) 



where the components Xq and z are the probability 
amplitudes associated with the basis vectors Eqs. (Cla) 
and (Clc) in the g and / sets, respectively. Hence in this 
basis the couplings g <-> e <-> / provide independent path- 
ways of excitation. Each g state is connected through a 
single pathway to a single / state. 



APPENDIX D: STOKES FIELD MS TRANSFORMATION MATRICES FOR THE J = 1 <-> 2 <-> 3 LINKAGE 



The Stokes field MS transformation yields three eigenvalues Xk of the matrix SS^ composed from the Stokes field 
coupling matrix of Eq. (45) 



,'1-k 1 
Afe = z + w cot [ — - — 7r + — arctan v 



(Dl) 



for k = 1, 2, 3, where 



-V146004 cos 12(9 + 857454 cos 89 + 2234532 cos 4(9 + 1524810 
4 



2u 



(839 + 909 cos 40) 

73002 cos 126> + 428727 cos i 



11 17266 cos 46» + 762405 



w = 



22960u+ 1932Oucos40 
709 cos 40 + 923 



14490 cos 40 + 17220' 
The Stokes field MS transformation matrix A is given by 



(D2a) 
(D2b) 

(D2c) 
(D2d) 



i<f>S 



pi A) (Ai)/n(^(Ai) p 2 A) (A 1 )/n(' 4 )(A 1 ) p^ ) (A 1 )/n^)(A 1 ) p^ ) (A 1 )/n( A )(A 1 ) 

Pi A) (A 2 )/n^(A 2 ) p 2 A) (A 2 )/n^)(A 2 ) p^ ) (A 2 )/n^)(A 2 ) pi A) (A 2 )/n^(A 2 ) 

Pl A) (A 3 )/n^)(A 3 ) p 2 A) (A 3 )/n^)(A 3 ) P 3 A) (A 3 )/n^)(A 3 ) pi A) (A 3 )/n^)(A 3 ) 

^V4 A) 4 A) /n { f 4 A) M A) 4 A) M A) 



(D3) 



where the polynomials (x) and the normalization n^' (x) read 



(A), 



Pi {x 

(A), 
P2 \ X 

P3 i X 

(A), 

Pa \ x 
n {A \x 



o i(2Vs-20 s ) 



sin 0(147OOx 2 - 980(2 + cos 29)x + cos 49 + 56 cos 29 + 63) 



24 
v/15 



,i(4>s-<t>s) 



cos 0(294Ox 2 - (308 + 280 cos 29)x + 3 cos 40 + 12 cos 29 + 9), 



sin 61(28(1 + cos 29)x - cos 49 - 4 cos 29 - 3) 



1 



cos 6»(1 - cos 46») 



sJ\p[ A \x)\ + |p 2 A) (x)| 2 + \p[ A \x)\ + |pf } (x) 



(D4a) 

(D4b) 

(D4c) 
(D4d) 

(D4e) 
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and the coefficients d\ A \ and the normalization are defined as 



d[ A) = -jWs- 2 * 3 ) cot 3 6 , 



(A) 



15 e *(*s-fcOcot 2 0. 



d{ A) = -VTEcote, 



d (A) = e^^s-V-s) ; 

n (A) = v 7 ! + 15 cot 2 + 15 cot 4 9 + cot 6 . 
Similarly, the other Stokes field MS transformation matrix is obtained as 



p{ B \x 1 )/n( B H\ 1 ) pi B \\ 1 )/n( B HX 1 ) p< B >(Ai)/n( B >(Ai) 
p[ B \x 2 )/n( B HX 2 ) pi B \x 2 )/n( B \X 2 ) p {B \X 2 )/n( B \X 2 ) 
p[ B \x 3 )/n( B HX 3 ) p {B \x 3 )/n( B HX 3 ) p 3 B \X 3 )/n^(X 3 ) 



(D4f) 
(D4g) 
(D4h) 
(D4i) 
(D4j) 



(D5) 



where the polynomials p\ D> (x) and the normalization n^ B \x) read 



p[ B \x) = e l ^ s -^ s) p[ A \x)/sm6 , 
p {B \x) 
p {B \x) 

n^(x) 



/fi 

= ^-sin26»(lO5a;-7cos20-8), 
= le^s-V-s) sin 2 20, 

= \I\p[ B \x)\ 2 + \A B) (x)\ 2 + \p ( 3 3 \x) 



(D6a) 
(D6b) 
(D6c) 

(D6d) 



(1) 


sin^e^ 


o — 


cosx 


(2) ... 


cos x e% ^ 


o — 


-sinx 



(1 2) 

The vectors x ' characterizing the dark states of Eq. (12) are obtained by finding the eigenvectors of the Hermitian 
matrix Eq. (14), which is obtained by inserting Eqs. (44), (47) and (D5) into Eq. (14). The two eigenvectors are given 

by 



(D7a) 
(D7b) 

(D8a) 
(D8b) 

(D8c) 
(D8d) 



where 



X 

v! 



1 2\u!\ 
- arctan — — 

2 v' 

argu', 
7 



s-i/>s) 



60 



sin 26(-8 + 7 cos 29 cos 2?y) 



+e 



i{<t>p-ip P ) 



sin 2rj 



7 



24 



343 7 



_ + _ _|_ _!_ e 2i(<^s-V's+V'P-0p) 



360 40 



sin 2 29 



49 



= — cos(0s — tps + ipp — 4>p) s i n 2»? sin 49 



301 203 



36 



90 



cos 29 I cos 2i] 



49 



cos 29. 



APPENDIX E: STOKES FIELD MS TRANSFORMATION MATRICES FOR THE J = 1 <-> 2 <-> 1 LINKAGE 



The Stokes field MS transformation matrix A is given by 



A = e^ s 



p[ A \x 1 )/n^\X 1 



Pi 



P2 



(A) (A 1 )/n^)(A 1 ) 



< A) (A 2 )/n^>(A 2 ) P {A) (X 2 )/n( A HX 2 ) 



(El) 



20 



where the polynomials p\ A \x) and the normalization n^ A \x) read 



p[ A \x) = -1-5 sin 2 (9 + 50a;, 
p {A \x) = sinflcosfle-W*-**), 



n^(x) 



p[ A \x) + p (A \x) 



Similarly, the other Stokes field MS transformation matrix is obtained as 



B 



B a 
B b 



where 



and 



sgn(sin 49 sin 9) 

sgn(cos 9) 



B h 



p[ B \\ 1 )/n^{\ 1 ) pf)(A 1 )/n( s )(A 1 ) pi B \\ 1 )/n^(X 1 ) 
p (B) (X 2 )/n^(X 2 ) p {B) {X 2 )/n( B )(X 2 ) p {B) {X 2 )/n^(X 2 ) 



,(B) ; (B) ,(B) I (B) JB) I (B) 

d\ '/n d d 2 '/n d d\ /n d 



(E2) 
(E3) 

(E4) 



(E5) 



(E6) 



(E7) 



( B) ( B} ( B) ( B) 

The polynomials p\ (x) with the normalization n d (x) ; the coefficients d\ with the normalization n d read 



p[ B \x) 
p {B \x) 
pf\x) 
n^(x) 



_ e -i(4>s-ips) „,.. 2 
^ ■ An 

sin 49 , 

4 



cos 6»(7sin 2 6» + cos^6l-50x) ; 



-i(i>s-4>s) 



sin 9(7 cos 2 9 + sin 2 : 9 -50a;). 



(B) 



4 B) - 



p[ A \x) + p {A \x) + pi A) (x) 



e^-V-s) sin 2 ^ 
sin cos , 

e ^(V- S -0 S ) cos 2 6/j 

Vl + sin 2 29 . 



(E8) 
(E9) 
(E10) 

(EH) 

(E12) 
(E13) 
(E14) 
(E15) 
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